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At the time of recombination, baryons and photons decoupled and the sound speed in the baryonic 
fluid dropped from relativistic, ~ c/%/3, to the thermal velocities of the hydrogen atoms, ~ 2x 10~^c. 
This is less than the relative velocities of baryons and dark matter computed via linear perturbation 
theory, so we infer that there are supersonic coherent flows of the baryons relative to the underlying 
potential wells created by the dark matter. As a result, the advection of small-scale perturbations 
(near the baryonic Jeans scale) by large-scale velocity flows is important for the formation of the first 
structures. This effect involves a quadratic term in the cosmological perturbation theory equations 
and hence has not been included in studies based on linear perturbation theory. We show that the 
relative motion suppresses the abundance of the first bound objects, even if one only investigates 
dark matter haloes, and leads to qualitative changes in their spatial distribution, such as introducing 
scale-dependent bias and stochasticity. We further discuss the possible observable implications of 
this effect for high-redshift galaxy clustering and reionization. 

PACS numbers: 98.65.Dx, 98.80.Es 



I. INTRODUCTION 

The early Universe was extremely homogeneous and 
isotropic, with small adiabatic density perturbations 
likely seeded during an epoch of inflation [i|-y| • The sub- 
sequent evolution of the Universe is well-described by a 
model containing baryons, cold dark matter (CDM), and 
a cosmological constant (A). This inflationary ACDM 
paradigm, with only six parameters, is simultaneously 
consistent with a wide range of cosmological observables 

One of the key features of the ACDM scenario is the hi- 
erarchical formation of structure: since the matter power 
spectrum A^ (k) is an increasing function of wavenumber 
k, the smallest perturbations collapse first, followed by 
their assembly into larger and larger structures. The for- 
mation of the first structures has recently become a major 
research area: aside from the intrinsic interest in under- 
standing the first galaxies, these objects are believed to 
be responsible at least partially for the reionization of the 
intergalactic medium [6|, |7| , and they are sensitive to the 
small-scale power spectrum of the dark matter, which is 
a powerful probe of dark matter microphysics. 

The evolution of density perturbations in the early 
Universe is generally described using linear perturbation 
theory, which treats overdensities and velocity fields as 
small quantities and hence neglects second order terms. 
Several previous works have extended this theory down to 
the post-recombination baryon Jeans scale [^, Q . Inter- 
est in direct observations of the high-redshift Universe via 
absorption in the redshifted 21 cm line jlO] has motivated 
more detailed investigation of the clustering of baryons 
during the epoch between recombination and reionization 
[ll|,|l2|, including the entropy and ionization fiuctuations 



in the baryons [13l - ll5l | . A deficiency of linear perturba- 
tion theory is that it does not describe the collapse of 
perturbations to form bound haloes, although analytical 
models such as the Press-Schechter formalism [16l.[l7l| are 
often used to estimate the halo mass function and clus- 
tering. In order to go beyond linear perturbation theory, 
one may use spherically symmetric Lagrangian hydrody- 
namic models Ts] or high-resolution 3D hydrodynamic 
simulations to follow the infall of baryons into the first 
haloes [I^ [23| ■ Since it is inherently nonperturbative this 
approach can, with incorporation of appropriate chem- 
istry and cooling processes, even be followed all the way 
to the formation of the first stars I2ll423ll . 



The principal purpose of this paper is to point out a 
new nonlinear effect in the growth of small-scale density 
perturbations that is active even at z ^ 1000. The idea is 
that prior to recombination, the baryons are tightly cou- 
pled to t he p hotons resulting in a standing acoustic wave 
pattern [2J|. Modern linear perturbation theory treat- 
ments including the CDM [2^ [2g| show a consequent 
relative velocity of the baryons and CDM since the latter 
does not suffer Thomson scattering and merely follows 
geodesies of the cosmic spacetime. At the time of recom- 
bination, the root-mean-square (RMS) relative velocity 
is 30 km s~^, and this is coherent over a scale of several 
Mpc comoving (the Silk damping scale [23] )• When the 
baryons recombine and are no longer tied to the photons, 
their sound speed drops to ~ 6 km s~^, and hence there 
is a highly supersonic relative velocity between baryons 
and CDM. This means that near the baryonic Jeans scale, 
perturbations in the baryons and CDM are advected rel- 
ative to each other in less than a Hubble time (and hence 
less than their growth time). This effect is investigated 
herein, and we find that it both suppresses the growth of 



small-scale structure, and leads to qualitatively new ef- 
fects in the clustering of the first bound baryonic objects. 

The suppression effect does not appear to have been 
present in previous analyses. Since it results from the 
coupling of large-scale and small-scale modes, it is non- 
linear and hence not present in linear perturbation the- 
ory. Since the large-scale modes involved are associated 
with the acoustic oscillations of the photon-baryon fluid, 
they are not properly modeled by hydrodynamic simula- 
tions whose box size is smaller than the acoustic horizon 
(~ 140 Mpc comoving). 

Understanding of the physical processes that deter- 
mine the collapse of the first dark matter halos and sub- 
sequent accumulation of baryonic matter in those halos 
is of paramount importance for interpretation of future 
data on reionization, high-redshift galaxies, and possibly 
dark matter substructure. In the present paper we intro- 
duce the formalism and focus on the key features of the 
relative velocity effect, leaving a detailed study of various 
applications for future work. 

This paper is organized as follows. In Section |TT] we in- 
troduce the effect and calculate the Mach number of the 
relative motion of dark matter and baryonic fluids at the 
time of recombination. We show that accounting for the 
relative motion leads to a suppression of the matter power 
spectrum near the baryon Jeans scale. In Section [IIII we 
compute the abundance and clustering properties of the 
first haloes taking account of the relative motion. The 
treatment is simple (it uses the linear Gaussian random 
field model for the large-scale density and velocity per- 
turbations in cells of size a few Mpc, and then uses an 
analytic model to compute the density of small haloes 
in each cell), but we believe it should capture the qual- 
itative result of the relative motion effect. We briefly 
summarize our results and outline possible future work 
in Section HVl 

The numerical results and plots shown in this pa- 
per assume a cosmology with present-day baryon den- 
sity fib.o = 0.044, CDM density fico = 0.226, dark 
energy density JIa^o = 0.73, Hubble constant Hq = 71 
kmsMpc"^, and adiabatic primordial perturbations of 
variance A^ = 2.42 x 10"^. 



II. GROWTH OF SMALL-SCALE STRUCTURE 

INCLUDING RELATIVE VELOCITY OF 

BARYONS AND CDM 

Before recombination, baryons are tightly coupled to 
photons via Thomson scattering and the sound speed 
is that appropriate for a radiation-dominated plasma, 
^ c/-\/3. Perturbations in the CDM component can 
grow, however, because the CDM experiences no drag 
against the radiation. As the universe expands and cools 
electrons recombine with protons and the universe be- 
comes transparent [28, 29| . This also leads to a kinematic 
decoupling of the baryons from the radiation, so that the 
baryons can fall into the potential wells created by the 



CDM. The effective redshift of decoupling is ^doc ~ 1020, 
which is slightly later than the surface of last scattering 
for microwave background photons because the baryons 
have lower inertia than the photons during this epoch 

M. 



A. Basic setup 

In the post-recombination gas, the baryonic sound 
speed is 






(1) 



where 7 = 5/3 for an ideal monatomic gas, /i — 1.22 
is the mean molecular weight including a helium mass 
fraction of 0.24, ttih is the mass of the hydrogen atom, 
and Tb is the kinetic temperature of the baryons. Here T^ 
is determined by a competition between adiabatic cooling 
and Compton heating from the CMB; we obtain it using 
the Recfast code [3l|, ^^ and parametrize it as: 



Tb(a) 



7cMB,0 



a/fli 



1 + (a2/a)3/2 



(2) 



with ai = 1/119, 02 = 1/115, and Tcmb.o = 2.726 K. 

While the baryonic velocity drops precipitously during 
recombination dark matter velocity remains unaffected 
and after recombination dark matter motion with respect 
to baryons become significant. The relative velocity can 
be written as: 



Vbc(k) = -[0b(k)-0c(k)], 



where k is a unit vector in the direction of k, and 
a^^V • V is the velocity divergence. 
The variance of this relative velocity is 



(3) 



(«bc(x)> 



dk 
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%ik)-e,{k) 
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^vbc(^)j 
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where A?(fc) — 2.42 x 10 ^ is the initial curvature per- 
turbation variance per Infc [33|. Integration of Eq. ^ 
at the time of recombination (z^ec — 1020) shows that 
dark matter moves relative to the baryons with rms ve- 
locities ~ 30kms~^ corresponding to a Mach number of 
M = Vhc/cs ~ 5. This supersonic relative motion allows 
baryons to advect out of the dark matter potential wells 
and significantly suppresses the growth of structure at 
wave numbers higher than 
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FIG. 1: The coherence scale of Vbc is determined by the 
range of scales over which Ajjjj.(fc) is non-zero. Here we plot 
Aji-,;.(fc), the variance of the relative velocity perturbation per 
In k, as a function of wavenumber k. The power spectrum 
drops rapidly at A: > 0.5Mpc~^, indicating that the relative 
velocity is coherent over scales smaller than a few Mpc co- 
moving. 



The relative contributions to Vhc from different scales 
are shown in Figure [TJ One can clearly see that there is 
no contribution to the relative velocity from the largest 
scales, which were outside the sound horizon at the time 
of decoupling, and that the dominant contribution arises 
from the acoustic oscillation regime, which have typi- 
cal velocities of a few times cA^ and suffer no Hubble 
damping at early epochs when pj ^ p^ [25|. At the 
smallest scales, the acoustic oscillations in the baryons 
are damped by photon diffusion, and the CDM veloci- 
ties are suppressed by Hubble drag during the radiation 
era. Thus we expect that Vbc contains contributions from 
scales ranging from the Silk damping length up to the 
sound horizon. This leads us to the conclusion that there 
is a separation of scales: the scales at which the first bary- 
onic objects will form (^ 10 kpc) is much smaller than 
the coherence length of the relative velocity field associ- 
ated with acoustic oscillations (few Mpc). This will be 
critical for our use of moving-background perturbation 
theory to follow early structure formation. 



equation for the gravitational potential (e.g. [3J]): 

dS 

-^ + a-Vc-V5c = -a-i(l + (5c)V-Vc, 

— -+a-i Vc-VVc = i?Ve, 

ot a 

-^+a-Vb-V5b = -a-i(l + (5b)V-Vb, 
ot 

^+a-(vb.V)vb = -^-i7vb 
ot a 



-a ^c^V(5b,and 



a~^V^$ = 47rGpni<5n 



(6) 



Here $ is total gravitational potential and subscripts c, 
b, and m stand for dark matter, baryons and total matter 
respectively. 

A more complete treatment would also follow the 
baryon entropy [l3| and ionization fraction |14| pertur- 
bations. We have not done this here, but we note that 
the moving background perturbation theory approach de- 
scribed here could be extended to accommodate these 
additional variables. 

The standard way to do large scale structure pertur- 
bation theory is to Taylor-expand in powers of the pri- 
mordial perturbations, e.g. Sc = Sc + Sc + Sc + ■■■■ 
One may then use the linear terms in the above equa- 
tion to describe the behavior at order n; for n > 2, the 
quadratic terms in Eq. ^ may be treated as a source for 
the order-n perturbation, written hierarchically in terms 
of orders < n ■35'-38]. This approach can even be ex- 
tended to include both baryons and CDM [33|. In our 
case, this is not appropriate: since there are relative bulk 
flow velocities between the baryons and CDM with Mach 
numbers of order 10, it follows that for perturbations at 
the baryonic Jeans scale the baryon and CDM compo- 
nents will be advected relative to each other by up to 
several perturbation wavelengths. Whether the standard 
perturbation series will converge in this case is an open 
question; even if it does, we expect that many orders in 
perturbation theory would be required. Therefore, we 
desire an alternative method to follow the growth of the 
smallest structures. 



Moving-background perturbation theory 
(MBPT) 



B. Fluid equations 

After recombination the small-scale inhomogeneities in 
the photons and neutrinos are rapidly washed out by free- 
streaming, and the dark energy is not yet dynamically 
significant. Also on small scales we can ignore the general 
relativistic (higher order in aH/k) terms. Thus we can 
write the evolution equations as the pressureless Navier- 
Stokes equations for the CDM, the Navier-Stokes equa- 
tions including pressure for the baryons, and the Poisson 



Our preferred method of following the earliest struc- 
tures is to do a perturbation analysis on a background 
where the baryons move relative to the CDM. The idea 
is that in the absence of density perturbations, but in the 
presence of a bulk relative velocity, there exists an exact 
solution to Eq. dH): 



(7) 



Vc(x,i) 


= <'\t), 


Vb(x,i) 


= ^^^yit), and 


$ 


= 6c^Sb = 0, 



,(bg) 



,(bg) 



where Vc and vj^ ' are constant with position and have 
temporal dependence ex l/a(i). Without loss of general- 
ity, one may boost to a different reference frame and set 
one of these (e.g. v^ ^') equal to zero. 

Since the relative velocity of the baryons and CDM is 
coherent over scales of several comoving Mpc, whereas 
the scales of direct interest for us are at a few baryonic 
Jeans lengths (~ 10 comoving kpc), the above moving 
background is an appropriate zeroeth-order solution in 
any small (~ 1 Mpc) region of the Universe. Thus we can 
imagine the Universe as composed of many individual 
patches, each of which has a different relative velocity 



,(bg) 



Small fluctuations on this background grow due 
to gravitational instability; their early stages of collapse 
can be modeled using linear perturbation theory around 
Eq. ([7]) using the local value of v^^^' . 

Perturbing around Eq. ([7]), and writing the perturba- 
tion variables Ub c , 



Vb(x,i)=v[,'^s'(<) + Ub(x,i) 



(8) 



and similarly for Uc, we may transform Eq. ^ into a 
system of equations involving Ub,c instead of Vb.c- Work- 
ing only to first order in the new perturbation variables 
{(5c, Uc, (5b, Ub, $} we find: 



dt + ""'"^''^ ■ ^^^ = 


= -a"^V • Uc, 


^^";+a-i(v^).V)uc . 


= ^Uc, 

a 


^+a-v^).V4 = 


= -a^^V-Ub, 


^;;+a-(v^).V)ub^ 


= -ffUb 

a 




—ac;gV (5b, and 



a^2V^$ = 47rG'p,n'5n 



It is convenient to transform these equations into 
Fourier space and use the last equation to eliminate $. 
We may also re-write the velocity equations in terms of 
the divergence 6, with Ui(k) — — mfc~^k0i(k) (i=:b or c), 
since under the approximation of barotropic flow of the 
baryons the vorticity remains zero until the development 
of structure formation shocks. We may also work in the 
bulk baryon frame, i.e. we may set 



v[,"s) = and v('^g) = -v\^f{t). (10) 

This reduces our system of equations to 



9(5, 



- L^i^s') .IrS -0 

„, — Vj^ KOc (7c 

at a 



-^rl"^> • k^c 



dt 

^r- = -t'b, and 
at 



iW 



[Qj, + rib^b) - 2He, 



50b 
dt 



3H^ 



(r!c<5c + r!b(5b) - 2H6b + 



c'k^ 



(11) 



Note that ric,b are evaluated at the appropriate redshift 
rather than taking on their present-day values. Our code 
evolves these equations, albeit with the scale factor a as 
the independent variable, which can be accomplished us- 
ing the replacement d/dt — aH d/da. It is important to 



,(bg) 



1/ a{t) when evolving 



note the time dependence v^^ 

Eq. dm). 

On large scales k <^ fcvbc ^ 40Mpc~^ used for galaxy 
clustering and even for Lyman-a forest studies, the Vbc 
terms in Eq. dill) are negligible. However, at A: > fcvbc, 
the advection terms become comparable to or larger than 
the Hubble expansion rate, and they must be taken into 
account. Note that this is true even if one's interest is 
only in the CDM perturbations, since the baryons con- 
tribute 17% of the energy density and hence their pertur- 
bations are important in Eq. (jlip . (As an extreme exam- 
ple, below the Jeans scale k > kj, the growth of structure 
in the CDM switches from the "standard" i5 ex o growth 
to a slower growth S ex aV25-24S2b/4-i/4 P^IJ^].) 



D. Small-scale transfer function and matter power 
spectrum 

The usual way to describe the small-scale distribu- 
tion of matter is to derive a transfer function T{k) that 
maps primordial to final potentials, and then to write 
the matter power spectrum Pm(fc), equal to the primor- 
dial power spectrum times |r(fc)p (times normalization 
factors [41|). We may solve the transfer functions includ- 
ing the relative velocity effect by solving the system of 
equations, Eq. PTj) . We evolve these from the redshift 
of recombination, where initial conditions are determined 
using Cmbfast [H], to z 
function, evaluated at z = 



= 40. The resulting transfer 
40, is clearly a function of 
(9) the local relative velocity 'v^^' and also of the angle i} 



(bg) 



be 

between the direction of the wave vector k and v^^ 

We may determine a local isotropically averaged power 
spectrum Pioc.m(fc; ^bc) by averaging over the direction 
of k, i.e. we write: 



,(fc;«bc) = /'cW 



1 



(5m (k; Vbc 



C(k) 



sm'dd'&. 



(^2) 

where C(k) is the primordial curvature perturbation and 
P^(fc) is its power spectrum. This power spectrum de- 
pends only on the magnitude of Ubc (we drop the super- 
script '^8-' to reduce clutter). In order to determine an 
overall effect on the small-scale matter power spectrum 
we need to average over a large number of coherence re- 
gions with different Wbc- The latter arises from linear 
perturbation theory on large scales and hence is well- 
described by a 3-dimensional Gaussian distribution with 
variance per axis 



vbc 



Vbc(x) 



(13) 
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EARLY HALOES 
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FIG. 2: Power spectrum of matter distribution in the first 
order CDM model (solid line) and with the Vhc effect included 
(dashed line) at the redshift of z; = 40. 



We may then average Eq. ()T2|) to obtain a globally aver- 
aged matter power spectrum. 



Intuitively, we expect the relative velocity effect to sup- 
press the small-scale power spectrum, since the moving 
baryons have pressure ^ Pb^^bc i^^ ^^^ CDM frame. This 
suppression is shown in Figure [2] where we plot A^(fc) = 
[A;^/(27r^)]Pm(fc) for the cases with and without the ef- 
fect of relative velocity. The power is most strongly sup- 
pressed around the Jeans scale fcj = aH/cs ^ 200 Mpc~^, 
where a difference of ~ 15% is computed. 

The effect of Vf,c is not limited to the suppression of 
power on small scales, but rather has an important im- 
plication for the distribution of the first bound struc- 
tures with respect to matter distribution as well as for 
the number densities of the first halos. To study these 
effects we ran a set of simulations in which the large- 
scale density and velocity fields were generated accord- 
ing to linear perturbation theory. We then used analyt- 
ical (Press-Schechter) arguments to predict the number 
of haloes formed in each cell of our cosmological box. 
This hybrid approach is computationally feasible on a 
single desktop computer since it does not have to numer- 
ically follow the small-scale modes, and should capture 
the rough magnitude of the effect. However, ultimately 
a simulation that follows the full nonlinear evolution of 
the small-scale modes will be required. The key reason 
for using approximate methods in the present study, as 
opposed to a full hydrodynamic numerical simulation, is 
our desire to introduce the concept of relative velocity ef- 
fect in the simplest and most intuitive way while allowing 
more detailed study to be performed by other research 
groups in an unbiased manner. 



We now investigate the formation of the first bary- 
onic objects, taking account of the relative velocity ef- 
fect. This is a difficult problem, which we only partially 
solve in this paper: one has acoustic oscillations in the 
photon-baryon plasma that travel ^ 140 Mpc, and si- 
multaneously one must resolve the baryon Jeans scale. 
We provide a computation based on the formalism de- 
scribed above: we generate a realization of a Gaussian 
random primordial perturbation on a 3D grid, and then 
to each cell we assign an overdensity 6\ (where the "1" 
referes to long-wavelength modes) using periodic bound- 
ary condition and a relative velocity Vbc derived from 
the linear density field. Initial values of 6i are obtained 
using the linear perturbation theory, as there is no sig- 
nificant difference between the theory with and without 
relative motion effect before the time of recombination 
when the values of Si are formed. Then, within each 
cell, we use the peak-background split to compute the 
number density of haloes. The new twist is that the 
small-scale power spectrum is modulated by the large- 
scale Vhc- (In some ways, this is similar to the modifica- 
tion of the peak-background split used for local /NL-type 
non-Gaussianity studies |43l. |44| . except that in our case 
the modulation of the small-scale power spectrum is a re- 
sult of the advection process and arises even in standard 
ACDM cosmology with Gaussian adiabatic initial condi- 
tions.) This of course depends on an analytic model for 
the halo mass function; we have used the Press-Schechter 
model [ia, [13] • The validity of Press-Schechter for any 
precise calculation is dubious - particularly since it is 
being applied here with an anisotropic local power spec- 
trum - but we expect that the qualitative results (a scale- 
dependent enhancement in the bias and stochasticity at 
large scales, with acoustic oscillations in each) would still 
arise in a more accurate treatment. 



A. Peak-background split 

The collapse of the first halos can be conveniently 
treated in the framework of the peak-background split 
45[, in which the growth of small-scale inho- 



formalism 

mogeneities depends on the large scale overdensity. One 
can split the density field into a long- wavelength piece 5i 
and a short- wavelength piece Ss'. 



p(x) = p[l + Ji(x) + 5,(x)]. 



(14) 



In any region, the number density of haloes of any given 
type generally depends on the large-scale overdensity (5i, 
and on the statistics of the small-scale perturbations 5s 
(in particular, their local power spectrum). In the usual 
case where the small and large-scale perturbations are 
independent, the number density becomes purely a func- 
tion of the large-scale overdensity plus a stochastic com- 



ponent e with (e(x)) = 0; Taylor-expanding in 6\ gives 



The bias is then 



n[l + 6o(5;(x)]+e(x). 



_i dn 



(15) 



(16) 



This argument leads to a generically scale-independent 
bias at sufficiently large scales (with the addition of a 
Poisson or halo-shot-noise term |46ll47|). 

When the relative motion of dark matter and baryons 
is introduced the growth of small scale overdensities be- 
comes dependent on the local value of the relative veloc- 
ity. Equation (fT5|) then generalizes to 

n(x)=n[(5i(x),Wbc(x)]+e(x). (17) 

At Mach numbers of order 10, it is not clear whether we 
can Taylor-expand in Wbc- Therefore, our strategy will be 
to re-compute n(x) in each cell, using the Press-Schechter 
conditional mass function, i.e. the number of haloes per 
unit comoving volume per InM: 



A^(M|,^i,«bc) 



P Ss. 



^/2^ <J^ 



da 



dM 



x(l + <5j)exp[-^^2^^^], (18) 



where ^sc is critical overdensity of spherical collapse, M is 
the halo mass and N{M\6i,Vhc) has units of Mpc~^.The 
factor 1 4- (5i is the conversion from the Lagrangian vol- 
ume element (in which the Press-Schechter formalism is 
native) to the Eulerian volume element. Here a^ is the 
variance of the density field smoothed with the top-hat 
window function. 



a2(A/,Vbc)= / A^(fc,Vbc)|VF(fc,i?)|- 



, dk 



(19) 



where for A^(fc,Vbc) we use the isotropically averaged 
local matter power spectrum. In principle, one should 
follow here the formation of haloes in a statistically 
anisotropic density field. This will ultimately require a 
hydrodynamic (or at least A^-body) simulation to achieve 
results that can be used for detailed analysis. However, 
for the moment we use the Press-Schechter formalism; 
the top-hat window function in Fourier space can be writ- 
ten as W{k,R) — 3ji{kR)/{kR), where the smoothing 
scale R is determined by the halo mass M = ^nR^. 

In our case - unlike the usual case of purely Gaussian 
density perturbations - a'^{M) and hence da/dM are ex- 
plicit functions of relative velocity and hence will change 
from place to place. 



B. Simulation parameters 

Our fiducial box size is 1365"^ Mpc'^. The box is di- 
vided into smaller boxes each of the size of coherence 



length for the relative velocity field, A = 3 Mpc. Initial 
(i.e. at Zdcc) density and velocity distributions are gener- 
ated using the Cmbfast power spectrum computations 
[42| and smoothed with a Gaussian window function with 
scale length fcgcaie = 7r/(-\/3 A). The smoothing is neces- 
sary to avoid aliasing and spurious effects from the finite 
resolution of the simulation. For each small box we gen- 
erate initial values of Vbc and large scale overdensity Si at 
the time of recombination using Cmbfast. Halo number 
densities in each cell are then inferred from Eq. (J18p . 

In Figure [3l we show an example output from this pro- 
cedure. The top panel shows the matter density contrast 
and the bottom two panels the halo density contrast for 
Mhaio = 10*^ Mq without (Middle panel) and with (Bot- 
tom panel) the relative motion effect at the redshift of 
z = 40. Note that the structures in matter and halo 
overdensities, while correlated, are not identical. Com- 
parison of the halo density contrasts for two different 
cases clearly shows the importance of relative velocity 
effect on the formation of first bound objects. 



C. Halo abundance 

To illustrate the effect of relative velocity on the abun- 
dance of small haloes, we calculate number densities of 
collapsed halos with and without relative velocity. The 
decrease in number density is quantified by 



N^ 



No 



^N 



No 



(20) 



where A'o and Ny^c are average number densities of halos 
without the effect of Vbc and with it. The comparison 
of the two cases shows that the number density of haloes 
is suppressed by more than 60% at the mass scale of 
M ~ 1O^M0, as can be seen in Figure ID Note that the 
strongest suppression occurs for halo masses of 10^-^ M©, 
corresponding to top-hat scales of 20 kpc comoving, i.e. 
near fc~bc- ^^ emphasize that the results provided in the 
figure are based on the Press-Schecter formalism and are 
a good qualitative guide, but should not be interpreted 
quantitatively. 



D. Bias, stochasticity and the large scale 
distribution of early haloes 

The introduction of relative motions modifies the cor- 
relation between the first halos and the matter distribu- 
tions rendering bias parameter scale dependent. Because 
of the non-linear terms in the evolution equations dark 
matter and baryonic matter evolve out of phase and the 
growth of the overdense regions become dependent on 
both (5i and Ubc- 

To quantify this effect we calculated halo overdensity 
using number densities of halos in each of the small boxes 
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FIG. 3: Top panel: The matter density contrast 5m on a 2D 
slice of the 3D simulation box. The halo density contrast 5„ 
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panel) and with Vjc 7^ (Bottom panel). All panels are at 
z = 40. 
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FIG. 4: Top panel: The number density of dark matter haloes 
produced in our simulation box without the effect of relative 
velocity (solid line) and with the effect (dashed line). Bot- 
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haloes as a function of the halo mass. Number densities in 
our simulation correspond to the redshift oi z = 40. 



from our simulation: 



5n{M,x) 



N{M,x)- N{M) 
N{M) ' 



(21) 



Next, we calculate power spectra of halos of various 
masses: 

(27r)3fc(k-k')Phh(fc|M) = {6„{M,k)6n{M,-k'r), (22) 

where du is the Dirac (5-function. 

The difference between this case and the case neglect- 
ing Wbc can be illustrated by defining a bias correction 
parameter Ab{k): 



Phh(fc) - b, 



1 



Ab{k) 



J^^^ram\'^ ) ; 



(23) 



where 60 is a Gaussian scale independent bias, which in 
the Press-Schechtcr formalism is given by: 



1 



+ 1. 



(24) 



Using these results along with the matter power spec- 
trum we can obtain the scale-dependent component of 
the bias parameter Ab which is plotted in the top panel 
of Figure [S] for various halo masses. The plot shows that 
for halos with mass M ~ W^-IO^Mq there is a signif- 
icant increase of the bias. The effect of Vbc becomes 
less important for heavier and lighter masses which can 
be expected from the analysis of power suppression in 
Fig. [2l This is principally a consequence of the fact that 
for very massive haloes the baryons advect through a dis- 
tance that is only a small fraction of the halo scale R, and 
hence this advection does not affect the formation of the 
halo; whereas for the lowest-mass haloes, whose scale R 
is smaller than the baryon Jeans length, the baryons can 
be treated as homogeneous irrespective of their velocity. 

To further understand the importance of Wbc we calcu- 
late the stochasticity of the halos relative to the matter. 
In the bottom panel of Figure [S] we plot the stochasticity 
X as a function of wave number k for various halo masses. 
The stochasticity is defined as: 



X 



PLik) 



Pbh{k)PnUk)'^ 



(25) 



where the cross power spectrum Phm is defined via 
(27r)3,5D(k-k')Phm(fc) = (5h(M,k)5;;(k')). In the linear 
theory, without consideration of the v^c effect one would 
have X = 1 (modulo Poisson corrections as described 
above) . 

We checked the convergence of our results by run- 
ning the simulation with varying box sizes and vary- 
ing A. Specifically we tried runs with A = 4 Mpc, 
and found changes of less than 1% in the stochastic- 
ity and bias over the range 0.2 < k < lMpc~^ at 
-^haio — W^Mq, whereas using A = 6 Mpc produces 
change greater than 5% and distorts functional forms of 
both bias and stochastisity at A: > 0.1 Mpc~^. Similarly, 
increasing the box size to 2275^ Mpc^, whith fixed A did 
not produce observable change in x a-nd Ab, whereas de- 
creasing the box size to 1000^ Mpc^ changes our results 
by ^ 5% at fc < 0.1 Mpc^^. As a test, we repeated the 
analysis setting Vbc = to recover the 'standard' picture 
with a scale-independent halo bias and the stochastisity 
consistent with the linear theory prediction. Specifically, 
we found that at fc < 0.2 Mpc~^ the stochastisity is 
0.98 < X < 1- The small deviation from unity can be 
explained by the fact that mapping from the overden- 
sity Si to the number density of halos iV(Af |^i, Wbc) is not 
exactly linear even in Press-Schecter model. 

We also would like to mention that Figure [5] exhibits 
strong oscillations of Ab which correspond to the BAO 
in the matter power spectrum. This means that the sig- 
nal of the BAO in the halo power spectrum for these 
halo masses is very different from that of the dark mat- 
ter. To illustrate this point we plot the actual scaled 
halo power spectrum Phh{k)/bQ{Mh) (Figure [6]) for dif- 
ferent halo masses covering the range from Af/,, ~ 10** to 
Mh — 10^ that shows the behavior of the BAO signal. 
Although these are very low-mass haloes compared to 
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FIG. 5: The correction to the bias parameter A6 (top panel) 
and the stochasticity x = ^hm (bottom panel) for various 
halo masses ai z = 40. The solid curve corresponds to Mh = 
10"* M©; the thick-solid to Mh = 10^ M©; the dashed to Mu = 
IffM^; the dot-dashed to Mh = IO^'Mq; and the dotted to 
Mh = 10* Mq. In the first order CDM model A6 = and 
X = 1 on large scales. The enhancement of bias on small 
scales k > 0.3 Mpc"'^ is due to the nonlinear dependence of 
halo abundance on 5i. 



those probed by BAO surveys (M > lO-'^^M©), they are 
the seeds of present day galaxies, and their subsequent 
evolution might alter the BAO signal in the galaxy power 
spectrum at lower redshifts. As with other interesting 
applications of the relative velocity effect we relegate de- 
tailed analysis of the this problem to a future study. 



IV. CONCLUSIONS AND PROSPECTS 

We have shown that the relative velocity of baryonic 
and dark matter fluids plays an important role in the 
formation and evolution of small-scale structure of the 
early universe. In light of the increasing interest in reion- 
ization, high-redshift galaxy clustering, and dark matter 
substructure, it is imperative to understand the evolu- 
tion of small scales structure and all physical effects that 
contribute to this evolution. Here we discuss the possi- 
ble implications and next steps in exploring the relative 
velocity effect. 

Early galaxies may be observed in the next decade, 
either directly via the James Webb Space Telescope or 
indirectly through the near-infrared background and its 
fluctuations. Due to the relative velocity effect, the 
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and the dotted to Mh = 10* Mq. 



large-scale clustering of these galaxies should show scale- 
dependent bias and (if sufficient statistics are available 
to split the galaxy population into multiple samples) rel- 
ative stochasticity between different samples of galax- 
ies. For example, in our calculation at z — 15, haloes 
of mass 10* Mq show an increase of A6 — 0.73 over 
the no-velocity result of bo ~ 4.6 at fc = 0.02 Mpc""'^. 
Whether this effect will be detectable depends on the as- 
yet-unknown luminosity function of the highest-redshift 
galaxies, and whether the relative velocity effect can be 
separated from the scale-dependent bias produced by 
reionization (48ll49|. 

Since reionization is believed to be driven by the for- 
mation of the first haloes massive enough to contain gas 
and produce stellar photoionizing radiation [50] , the rela- 
tive velocity effect will delay reionization. However, this 
effect is degenerate with the mapping from halo mass 
to the number of massive stars formed and given the 
modest (factor of ~ 2) effects investigated here we do 
not expect that the effects can be disentangled from the 
reionization history alone. More interesting would be an 
investigation of the spatial structure of reionization and 
of the related high- redshift 21 cm signal, which has been 
investigated analytically and in simulations |5ll456l | . The 
scale-dependent bias and stochasticity we find here for 
early haloes may have a significant effect on the struc- 
ture (power spectrum, topology, and correlation with the 
matter density field) of early reionization bubbles. 

If the pre- reionization 21 cm signal jlO[ is ever ob- 
served, and cosmological information extracted, the rel- 
ative velocity effect will be very important: the smallest 
scale fluctuations in the baryons are modified at the tens 
of percents level. Indeed, since the 21 cm signal is non- 
linear in the baryon density (in the limit where the hy- 
drogen spin temperature is closely coupled to the CMB 
temperature, the signal is proportional to n? times the 
temperature-dependent collision cross section [57[), it is 
likely that even the large-scale fluctuations would be af- 
fected because regions with increased small-scale baryon 



power spectra will show more absorption. The locally 
anisotropic nature of the small-scale baryonic perturba- 
tions would also represent an issue for weak lensing of 
the 21 cm field |58l - l6l| and/or non-Gaussianity searches 
[62, |63|- A full analysis of the effect on the 21 cm power 
spectrum and non-Gaussian statistics is deferred to fu- 
ture work. 

Finally, the ACDM cosmology predicts that early dark 
matter haloes in the affected range of scales (mainly 
^ lO^'-lO^Mo) are assimilated into larger structures. 
Some of these early haloes may still be present today as 
dark matter substructure, which has attracted a great 
deal of interest since the subhalo mass function is in 
principle sensitive to the primordial small-scale CDM 
power spectrum and hence to possible deviations from 
"vanilla" CDM behavior (e.g. warm dark matter, or 
particles that are kinetically coupled to baryons at high 
redshift). Unfortunately, the overall power suppression 
effect we describe is probably not detectable via sub- 
structure since the transition from an initial CDM power 
spectrum through the formation and survival of substruc- 
ture is still not quantitatively understood (e.g. t64"'66j ) . 
However, the power suppression effect is modulated by 
the relative velocity field, which comes primarily from 
large-scale modes in the primordial density fluctuations 
(fc ~ O.lMpc"^) and hence can be reconstructed from 
large-scale structure surveys. Therefore it would be valu- 
able for future work to investigate whether the Ubc effect 
can be detected by differential measurements that com- 
pare the substructure abundance in strong lens systems 
where the lens haloes have similar mass but different re- 
constructed Vbc- It is also important to mention that 
the suppression of the formation of the flrst halos, which 
seed present day galaxies, and the decrease in the high-k 
power spectrum might help alleviating the known prob- 
lem of the over-abundance of substructure of dark matter 
halos in the ACDM model (the missing satellite prob- 
lem). This effect might also be important for predictions 
of the annihilation signal from dark matter particles. We 
relegate detailed investigation of these questions to the 
future study. 

In summary, we have shown that in the post- 
recombination Universe, there are bulk relative motions 
between the baryons and dark matter that are supersonic 
and are coherent over scales of several comoving mega- 
parsecs. The combined growth of small scale structure 
(between the baryon Jeans length Aj and A^Aj, where 
A4 is the Mach number of relative motion) is suppressed 
due to the baryons advecting out of the potential wells 
created by the dark matter. We find at lower redshifts 
(e.g. z ~ 40) a suppression of the power spectrum by 
^ 10% on scales of 50-500 Mpc^^ that is highly spatially 
variable. The suppression results in some reduction in the 
abundance of early haloes, but more importantly changes 
their spatial structure, leading to scale-dependent bias 
and stochasticity of the first haloes. These latter effects 
may be large for early low-mass, high-bias haloes, e.g. we 
find squared correlation coefficients x = ^hm ^^ small as 
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^ 0.2 at z = 40. Whether this unusual clustering pattern 
affects the spatial morphology of reionization depends on 
the importance of low-mass haloes and hence is unknown 
at this time, although we note that future 21 cm obser- 
vations combined with simulations that distribute their 
sources of ionizing radiation in different ways may shed 
some light on this issue. Farther in the future, the power 
suppression effect would certainly be significant for the 
interpretation of any pre-reionization 21 cm signal. In 
any case, our analysis highlights the importance of recon- 
sidering "standard" notions of structure formation (e.g. 
linear bias of haloes on large scales) as we enter new 
physical regimes at high redshift. 
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